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Abstract 

Two integrable differential-difference equations are derived from a (2-|-l)-dimensional modified Heisen¬ 
berg ferromagnetic equation and a resonant nonlinear Schroinger equation respectively. Multi-soliton 
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1 Introduction 

Recently, integrable discretizations of integrable equations have been of considerable and current interest 
in soliton theory. As Suris mentioned [l|, various approaches to the problem of integrable discretization 
are currently developed, among which the Hirota’s bilinear method is very powerful and effective. Discrete 
analogues of almost all interesting soliton equations, the KdV, the Toda chain, the sine-Gordon, etc., can 
be obtained by the Hirota method. The purpose of this paper is to consider integrable discrete analogues of 
two nonlinear Schrbdinger (NLS)-type system by Hirota method. 

The one-dimensional classical continuum Heisenberg models with different magnetic interactions have 
been settled as one of the interesting and attractive classes of nonlinear dynamical equations exhibiting the 
complete integrability on many occasions. As is well known, Heisenberg first proposed in 1928 the following 
discrete (isotropic) Heisenberg ferromagnetic (DHF) spin chain Q 

S„ = S„ X (S„+i + S„_i), (1) 

where Sn = (s", 82 , 83 ) G with |Sn| = 1 and the overdot represents the time derivative with respect to t. 
The DHF chain plays an important role in the theory of magnetism. 

A performance of the standard continuous limit procedure leads DHF model © to the integrable Heisen¬ 
berg ferromagnetic model 


St = S X (2) 

which is an important equation in condensed matter physics Q. NLS-type equations are extensively used 
to describe nonlinear water waves in fluids, ion-acoustic waves in plasmas, nonlinear envelope pulses in the 
fibers. It is known that HF is gauge equivalent to NLS equation and DHF is gauge equivalent to a kind of 
discrete NLS-like equation Q. 

Higher dimensional nonlinear evolution equations are proposed to describe certain nonlinear phenomena. 
Due to the dependence on the additional spatial variables in higher dimensional systems, richer solution 
structure might appear, such as dromions, lumps, breathers and loop solitons. 
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An extension of Eq. © is the (2 + l)-dimensional integrable modified HE system as follows, 


Ut + Uxy +UW = 0, 

(3a) 

Vt - Vxy -VW = 0, 

(3b) 

Wx + {uv)y = 0, 

(3c) 

which is associated with the (2 -|- l)-dimensional NLS equation 

i9r + gjT - 2g / (|gpg)r,d^ = 0, 

(4) 


where 77 is a spatial variable w. System j[3]) can also be used to model the biological pattern formation in 
reaction-diffusion process [1,[^ . In Ref. Q, system ([3|) was investigated through the prolongation structure 
and Lax representation. In Refs. 0,1^, integrable property of system ([3]) was studied through the Painleve 
analysis, and some localized coherent and periodic solutions were given by means of the multi-linear variable 
separation approach. Multi-soliton solutions of system (|3]) was derived in 0 by means of the Hirota bilinear 
method, and the double Wronskian solutions was given therein. Similar as the counterpart between DHF 
0 and continuous HE m, it is natural to consider discrete version of the (2 -|- I)-dimensional modified HE 
system ®. 

A new integrable version of the nonlinear Schrddinger (NLS) equation, called the resonant NLS (RNLS) 
equation. 


iUt 


a,„,2T- ^ 1 


■Uxx+ 4 |U|^U = /3p|UUxU, 


(5) 


was recently proposed ll| to describle low-dimensional gravity (the Jackiw-Teitelboim model) and response 
of a medium to the action of a quasimonochromatic wave. Here a is a nonlinear coefficient, and /3 denotes 
the strength of electrostriction pressure or diffraction. The term |C/|a;a,/|C/| on the righthand side of Eq. (O 
is so called the ’’quantum potential”. Moreover, Eq. ([5]) can model propagation of one-dimensional long 
magnetoacoustic waves in a cold collisionless plasma subject to a transverse magnetic field 0. Note that 
when /3 < I Eq. ([5]) is reduced to NLS 


( 6 ) 


-I- ^rr + Cr|$p5' = 0. 

When /3 > 1, it is not reducible to the usual NLS equation but to a reaction-diffusion (RD) system [H 

The Lax pair in the 2x2 matrix form of Eq. ((S]) was given in . Based on the RD version of RNLS, multi- 
soliton solutions were derived via the Backlund-Darboux transformations in [l^ and via the Hirota method 
in 0 respectively. In 141, integrable extension of Eq.® has been suggested and soliton solutions of the 


resonant NLS case under the reduction condition have been presented via the Hirota method. Additionally, 
the method of binary Bell polynomials was used to study two-soliton solution and integrable properties of 
Eq.® in [l6| . 

Since NLS equation has integrable discrete versions and RNLS is an intermediate equation between the 
focusing and defocusing NLS equation, it is natural to study the discrete analogue of RNLS. This paper is 
to investigate the integrable discretization of RNLS equation ® in the case of /3 > 1. 

The paper is structured as follows. In section 2, we present multi-soliton solutions and discrete version of 
(2 -|- 1)—dimensional modified HE equation. In section 3, we study RNLS equation, including multi-soliton 
solutions, semi-discrete analogue and dynamic properties of solutions. Finally, a short conclusion is given in 
section 4. 


2 (2+l)-dimensional modified HF equation 

2.1 iV-soliton solution to modified HF system 

Through the variable transformation 

w = 2{\n{F))^y, (7) 
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Eq.Q transforms into the following bilinear form 


{Dt+D,Dy)G,F = 0, 

( 8 a) 

iDt-D,Dy)H,F = 0, 

( 8 b) 

DlF,F + GH = 0. 

( 8 c) 


The Hirota bilinear differential operator is defined by 0 


m, n = 0 , 1 , 2 , • • • . 


/_a 

dx') \dy 



n 

a{x,y)b{x',y')\j:'=x,y'=y, 


One-, two- and three-soliton solutions of bilinear equation set (jS)) are given in 0 - Note that the bilinear 
equations (l 8 |) are in Schrbdinger type, here we present a compact form of multi-soliton solutions to the 
system dH) . The two-soliton is given as 


F = 1 -I- a(l, 1*) exp(7?i -I- ^i) -I- a(l, 2*) exp(77i -|- ^ 2 ) 

-h a(2,1*) exp(772 -H 6) + a(2,2*) exp{r]2 + 6) 

-I- a(l, 2,1*, 2*) exp(r 7 i -I- 772 -f Ci + ^ 2 ), 

G = exp(? 7 i) -H exp(r 72 ) -H a(l, 2,1*) exp(? 7 i + t ]2 + 6) 
+ a(l, 2, 2*) exp(? 7 i + 1 ^ 2 + 6), 

H = exp(^i) -H exp(^ 2 ) + a(l, 1*, 2*) exp(? 7 i + + ^ 2 ) 

-h a( 2 , 1 *, 2 *) exp (?72 -H Ci + 6 )- 

with 


r]i = kiX + ky + LOit + rfl, Wi =-k^k, 
— Pi^ T Qiy -f Ojt -t- 


and the coefficients are defined by 


— 2 (/ti /tj) , 

1 


^(^15^2;*** 7^n) — 


liihik)- 


l<Z<fc<n 

Generally, exact 7V-soliton solution to Eq. (|S]) is expressed in the following form 


(0 

F= ^ exp 

p^ 0 ,l 

(o) 

G = ^ exp 
1 ^= 0,1 


N 2N 2N 

PjVj + Pjij-N + pipjAij , 

j=l j=N+l l<i<j 

N 2N 2N 

^3^3-N + k'ih'jAij , 

3 = 1 j=N+l l<i<3 


(o) 

H= ^ exp 
A= 0,1 


N 2N 2N 

Xj^j-N + XiXjAij , 

i=l 3=N+i l<i<i 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 


(18) 

(19) 

( 20 ) 
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with 


rij = kjx + IjU + ojjt, u!j = —kjlj, j = 1,2, ■ ■ ■ , N, (21) 

~ Pj^ “1“ QjU “I" ^jt, rij = PjQj, J = 1; 2, • • • , N, (22) 

exp(Ay) =-2(fcj - fc^)^, i < j = 2,Z,-■ ■ ,N (23) 

exp(A,,jv+j) = +p-)2 ’ ,iV, (24) 

eyiY>{AN+i,N+j) = -‘^{Pi-Pjf, i < j = 2,S,-■ ■ ,N. (25) 


Here aj, "fj are both real parameters relating respectively to the amplitude and phase of the ith soliton. The 
sum 1 indicates the summation over all possible combinations of /x,; = 0,1 under the condition 

N N 

^ ^ ^ pN+j 1 (26) 

i=i 1=1 

and 1 ’ Sa=o i indicate the summation over all possible combinations of Xi = 0 , 1 , Vi = 0,1 under the 

condition 

N N 

^ ^ i^N+j + 1, (27) 

1=1 1=1 

N N 

E^1=EW-1- (28) 

1=1 1=1 

The proof of the Ai-soliton solution here is similar to that of the combined Schrodinger-mKdV equation in 
[isl and can be completed by induction. One can check the details therein. 


2.2 Integrable semi-discrete analogue of modified HF equation 

In this section, we construct the integrable discretization of HF equation ([5|) by using Hirota’s discretization 
method M- The Hirota bilinear difference operator exp((5Il„) is defined as (l7| . 


exp((5Zl„)a(n) • 6 (n) = exp 


Y ^ d 

* V dn dn' 


a(n)b{n') a{n + 5)b{n — S). 


The differential-difference HF system is obtained by discretizing the spacial part of the bilinear Eq. 


DxG*F 1 —(Gn+lFn — GnFn^i), (29) 

e 

DlF,F^^ (E„+iE„_i - Fl) , (30) 

where x = ne, n being integers and e a spatial-interval. Substituting (Oni) - (l5ni) into Eq.® results to 

DtGn»Fn + —Dy{Gn+l • Fn — Gn* Fn+i) = 0 , (31) 

DtHn •Fn — —Dy{F[n+l • Fn — Fin • Fn+l) = 0, (32) 

^ {Fn+lFn-1 - Fl) + GnHn = 0. (33) 

We demand that the discretized bilinear forms are invariant under the gauge transformation: 

Fn -1 Fn exp(gon), 

Gn G„exp(go?T-), 

Hn Hnexp{qon). 
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(34) 

(35) 

(36) 






Then we find a gauge invariant semi-discrete bilinear HF equations 

{eDt + 2Dy)Gn+i • Fn + (eZ?( — 2Dy)Gn • F^+i = 0 , 
{eDt — 2Dy)Hn+i ^Fn + {eDt + 2Dy)Hn • F'ri+i = 0, 


- {Fn+lFn-l - F^) + GnHn = 0 , 


or equivalently in compact form 


{eDt cosh ^ -I- 2 Dy sinh ^)G'„ = 0, 

{eDt cosh ^ - 2 Dy sinh ^)i/„ . = 0, 

2 (F„+iF„_i - Fl) + e^GnHr, = 0 . 


Let 






Un = 


Vn = 


r„ = In ■ 


F 


n+1 


p p 7 p 

n n n 

Then Eqs. (I37D - (I39P are transformed into ordinary nonlinear form 

2 r 

{Un+1 + Un)t + {Un+1 “ Un)rn,t + “ {Un+1 “ Un)y + (Un-l-l + Un)r, 

2 r 

{Vn+1 + Vn)t + ('Cn-l-1 “ -(u„+i - Vn)y + {Vn+1 + Vn)rn,: 


e L 


= 0 , 
= 0 , 


r(e^ 


- 1) -I- UnVn = 0. 


(37) 

(38) 

(39) 


(40) 

(41) 

(42) 

(43) 

(44) 

(45) 

(46) 


When we take the continuum limit e —)• 0, (jll9ll - (ll29() reduce to (l77])-(l78|) and 

2{\ii.F)^x + uv = Q. (47) 

Differentiating Ea. ()130p with respect to variable y, we get Eq. ([3c]). Thus we regard (I119ll - (|129|1 as a semi¬ 
discrete version of the HF system ([3]). In the following discussion we take the interval e = 1 for the sake of 
simplicity. 

Following the Hirota method, we expand and Fn in series with a small parameter 6 as 


= 1 + 5^f[^^ + 

+ • • • + + • • • , 

(48) 

Gn = 5G£i) + (53g£3) + .. 

.^^(2fc+l)^^2fe+l)^... , 

(49) 

Hn = 5h!;^^ + + ■ ■ 

.-4 j(2fe-ei)^^2fe+i) _^ 

(50) 


Substituting the expansion into the above bilinear Eqs. (EZI)-®, we find that there are only odd order 
terms of 6 in the first two equations while only even order terms appear in the third one. By the standard 
direct perturbation method, we obtain the one-soliton solution 


G„ = 7 iexp( 77 i), Hn ='y[exp{r][), Fn = 1 - exp{T]i + r][), 

where iji = pit + qiy + ln(/3i)n, ij'i = p[t -|- q[y -|- ln(/3()n, and /3i, /3( satisfy 

2 qi - Pi 


Pi = 


2 qi + Pi ’ 
2q[ + p[ 


^ ‘^q'l-p'i 

ai, 7 i and a'i, 7 ( are arbitrary constants. The two-soliton solution is presented as 

Gn = 5 [exp(? 7 i) -I- exp(r 72 )] -|- < 5 ^[ai 2 i exp( 77 i 4-772-1- p[) -f 0122 exp( 77 i 4-772-1- 773)], 
Hn = S[exp{T][) + exp( 772 )] -f S^[bi2i exp{p[ + P2 + Vi) + ^122 exp( 77 ( + P2 + V2)], 


Fn= 1 - 5 ' 


7i7l/3i/31 exp( 771-4770 717 ^^ 1 ^^ exp( 771 -4 77 ^) 727(/32/3i exp(772 -4 77 ^) 


2{PiP[ - 1)2 
I2I2P2P2 exp(772 -4 77^) 


mP2 - 1 )' 


2(/3i/3' - 1)2 

+ (5'‘xi2i2[exp(77i -4 772 -4 77( -4 772)] 


2(/32/3i - 1)2 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 
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where 


0121 = 


6121 = — 


X1212 = 


7i727i(/3i - 1^2^Pi 

'2(/32/3;-1)2(/3i/3J-1)2^ 

717(72 - P 2 ?Pl 


ai 22 — 


&122 — — 


717272 (/ 3 l - P 2 fP 2 

2 (/ 32 / 3 ^- 1 ) 2 (/ 3 i/ 3 ' - 1 ) 2 ’ 

727(7((/3( - P2fPl 


2(/3i/3; - l)2(/3i/3' - 1)2 ’ 2(/32/3^ - lY[P2P'i - 1)^ ’ 

{Pi - P 2 f{P'l - P 2 f-tl-f 2 l'll 2 PlP 2 P'lP 2 


4(/3i^; - l)2(/3i/3' - iy{P2P[ - 1)2(/32/3( - 1)2 • 

We can use the following compact expression for the above two-soliton solution, 


Fn = l + a(l, 1 *) 7 i 7 ( 7 ( exp(77i + r][) + a(l, 2^)717' exp(77i + 77^) 

+ a( 2 ,l*) 727 (exp (?72 + ? 7 () + a( 2 , 2 *) 7272 exp(r 72 +772) 

+ a(l, 2,1*, 2*)7 i 7(727( exp(r7i + 772 + r7( + 772 ), 

Gn = 71 exp( 77 i) + 72 exp(772) + a(l, 2 , l*)7i727( exp( 77 i + 772 + 77^) 
+ a(l, 2 , 2 *) 7 i 7272 exp( 77 i + 772 + 772), 

Hn = I'l exp( 77 () + 7^ exp( 77 ^) + a(l, 1*, 2*)7i7(7^ exp( 77 i +77^+ 77^) 
+ a( 2 , 1 *, 2*)727(72 exp(772 + 77 ^+ 77 ^), 


where the coefficients are defined as 




a{iG*) = - 


PiPj 

PiP' 


(ft/ 3 '- 1 ) 2 ’ 

f\2 


*(ft,r) = -2 


iPi-Pj) 

m 


( 57 ) 


(58) 

(59) 

(60) 

(61) 

(62) 

(63) 


and a{i,j, k*),a(i, j*,k*),a(i, 7, k*,1*) satisfy the operation rule (HB. In the same way, we can construct the 
three-soliton solution as that in [20j. The above expressions of the one- and two-soliton solutions suggest 
the exact Wsoliton solution of Eqs. (I37l) - (l3^ in the following form 


where 


ie) 

N 

2N 

2N 



Pn = y ] exp 

+ 

PjPj-N + 

y ^ PkPiAki 

5 

(64) 

M^0,1 

. 7=1 i 

^AT+l 

l<k<l 



( 0 ) 

N 

2N 

2N 



Gn = y ] exp 

+ 


y ] vkviAki 

7 

(65) 


i=i i 

^AT+l 

l<k<l 



( 0 ) 

N 

2N 

2N 



Hn = y ] exp 

y ^ + 


y ^ XkXiAki 

7 

(66) 


1=1 7 

^AT+l 

l<k<l 



Vj = Pjt + QjV + ln(/3ftn -h Jj, 

q 2gj - pj 

^ 2qj -1- pj 

J = l,2,--- 

,3V, 

(67) 

P'j = P'j^ + + MPj)n + I'v 

0 , 2g'+p' 

^ 2g'-p' 

J = l,2,--- 

,3V, 

(68) 

exp(Afe7) =fe</ = 2,3,--- ,iV 

PkPi 


(69) 


(R> _ /5''\2 

exp(^^+fc,^+ 7 ) = fc < / = 2,3, • • • , iV. 

PkPl 


(70) 

(71) 
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Here pj,qj, I 3 j,^j are all real parameters. The summations ii Si'=o i Si=o i satisfy the condition 

(| 26 )l - p^ respectively. 


2.3 Soliton propagation and interaction 

By virtue of (H 51 ) . (IM)) - (l 571 ) . interaction between the two solitons can be investigated. From Fig. 1 , head-on 
elastic interaction between the two solitons is found. With the time evolution, two solitons travel towards 
each other and then separate. After the interaction, the two solitons maintain their original amplitudes, 
velocities and shapes except for the phase shifts. To interpret the elastic interaction behavior of two solitons 
under the condition 71727172 ^ 0 , asymptotic analysis is carried out in the following. 

First we take the notation Xij = and let 

A = {| 77 i + v'i\ ,\m + ii2\ Av 2 + v'i\, \v2 + vA > \m - V2\, Wi - vA} ■ 


When 71727^72 0 , we have the following asymptotic behaviors as 


</> = 


GH 

yr 




| 7 i 7 iXii 
717 i 


|7i7ixii 
72 72 


172 72X22 

72 72 


172 72X22 


sech' 


2 / 7i7(xi2i2X22 


sech^ ( 


sech' 


2 / ^2+y?2+In|7272Xl212Xll^ I 


r]i + r][ fixed, 772 + ^2 —00, |?72 + ^72! = max A 

771 + 77^ fixed, 772 + 772 ^ +CX), I772 + ^72! = max A 

772 + 772 fixed, 771 + r][ —00, \rii +77^! = max A 

772 + 772 fixed, 771 + 77'^ ^ +CX), I771 + 77^1 = max A 


7172 


|7i72Xi2 

7172 


|7i 72X12 

72 7 i 


|727iX2i 

72 7 i 


|727iX2i 


f r?i-erj'-Hn|7i7'xi2| 


sech' 


2 ( m+’ri'2+MlW2Xl‘212X21 \ 


sech^ 

2 ( ^ 2 -eTjl-Hn| 727 ixi 2 l 2 Xr 2 ^ I 


sech' 


r]i + 772 fixed, 772 + Vi \V2 + Vi\ = max A 

771 -I- 772 fixed, 772 -I- 77'^ ^ - 1 - 00 ,1772 + ? 7 il = max A 

772 -I- 77^ fixed, 771 -I- 772 ->■ -00,1771 -I- 77'I = max A 

772 -I- 77'^ fixed, 771 -I- 772 -y - 1 - 00 ,1771 -I- 77^1 = max A 


( 72 ) 


71720122 

sech' 

I7172X12X22I 

71720121 

sech' 

I7172X11X21I 

7 i 72 *'212 

■sech' 

| 7 l 72X21X22! 

7 l 72*'112 

■sech' 

| 7 l 72X11X12 1 


»7l-i72-l-ln|7iXl272 '■X22’'| 
2 

m-'^2-l-ln|7iXll72”''x7i’'I 
2 

»7l-i72+ln|'7lx2l72“^Xi2^ 

2 

»7l-'^2+ln|7lxil72“^xr2^ 

2 


771 — 772 fixed, 77^ — 772 —)■ —00, | 77 ( — 772 1 = max A 

771 - 772 fixed, 77 ( — r]2 ^ -too, [772 + 77(1 = max A 

Vi ~ 772 fixed, 771 — 772 ^ —00,1771 -I- 77(1 = max A 

Vi ~ 772 fixed, 771 — 772 ^ -foo, I771 -I- 77(1 = max A 


Plots of two solitons interaction at different times are shown in Fig. [T] 

In addition, inelastic interaction between the two solitons can be found. When 7^ or 7' in dSll)-® is 
chosen to be zero, soliton fusion and fission may occur during the interaction. Specially, when 72 = 0 and 
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t= -30 


/=0 



Fig. 1: Two-soliton interaction for cj) = GH/F'^ by ((551) - (l501l with parameters 71 = 1, 72 = 1, 7 i 
— 1, Pi = 0.2, p2 = 0.3, qi = —0.7, q2 = 0.6, p'l = 0.3, P2 = 0.2, q[ = —0.8, 92 = 0-75 


-1, 72 = 


7i727i ^ Oj '"^6 have the following asymptotic analysis for (/), 


GH 


7 i 7 i« 


Vg ('/ 2 +i; 2 ) | ^27^6 + 7i727l^ai2ie''i ’’2 


{ vi+'n{+ri2+'n2 vi+'ni-V2-'n2 . ^ 

e 2 2 +X2i727ie 


' 717 I „„„-u 2 »?l+rjl+lnj7i7ixil 


sech^ 

7l7iXll| 

7271 1,2 »?2+rjl+ln|727lx2l 

7271X211 

7l72ai2l gggj^2 { r,i-i)2+ln|7lXll72~^X^i^| 


771 + r][ fixed, 772 + ?72 — t —00 

772 + 77^ fixed, 771 + 772 —>• —00 
771 — 772 fixed, r][ — ri2 ^ +00 


, I7172X11X21I 2 

Fig. 2 depicts one soliton split into two solitons with time evolution under the condition 72 = 0. 


3 Resonant nonlinear Schrodinger equation 

3.1 multi-soliton solutions of RNLS equation 

In [ 1 ^, only two-soliton solution of RNLS equation are given. In this section, we present N— soliton solution. 
We introduce the transformation 


U = exp 


u + iau + V — iav 


(73) 


where u = u{x,t) and v = v{x,t) are both real functions, and a is a real constant. Eq. ([5]) is transformed 
into 

2(7 — a)ut + (I — /3 — + 2ia)u^ + 2(1 — /3 + ia)uxx + 2(7 + a)vt 

+ (1 — /3 — — 2ia)Vx + 2(1 — /3 — ia)vxx + 2(1 — /3 + a^)uxVx + a exp(77 + u) = 0. (74) 

Since all the functions and parameters in Eq. (1741) are real, separating the real and imaginary parts of Eq. 
gives as 


2aut - (1 - /3 - a'^)ul - 2(1 - l3)uxx - ‘^avt - (1 - /9 - a^)vl 
- 2(1 - I5)vxx - 2(1 - /3 + a^)uxVx - a exp(77 + u) = 0, 
ut + aul + auxx + vt - avl - avxx = 0 . 


(75) 

(76) 


In what follows, we set /3 = 1 + a^. Adding and subtracting (1751) and (175]) yield, respectively, two similar 
equations, 

2(1 - a)[ut + aul + auxx] + 2(1 + a)[vt — avl “ aVxx] + = 0, (77) 

2(1 + a)[ut + aul + ^’^xx] + 2(1 — a)[vt — avl ~ o^Vxx] — = 0. (78) 

It is easy to check that provided Eos. (1771) - (1781) hold automatically, the following two equations are satisfied. 


a 


Ut + a{uxx + al) — — 0 ) 

4a 

Vt - a{vxx + vl) + = 0. 

4a 


Through the dependent variable transformation 


,9 . h 

u = in j, V = in j, 


Eqs. (|7^ - (|M1) are transformed into 

4a 


gj[iDt + aD^)g.f] - —(4a D^f.f + agh) = 0, 
0[(A - aDl)h,f] + j^{Aa^Dlf,f + agh) = 0, 


(79) 

(80) 

(81) 

(82) 

(83) 
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t= -30 


/=0 



Fig. 2: One soliton splits into two solitons for cj) = GH/F"^ by (1551) - (15(1)) with parameters 71 = 1, 72 = 1, 7 ^ = 
-1, 72 = 0, Pi = 0.2, p 2 = 0.3, qi = -0.7, 92 = 0.6, p[ = 0.3, P 2 = 0.2, q[ = -0.8, 92 = 0-75. 
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that can be decoupled into the following bilinear equations 


{Dt + aDl)g,f = 0, (84a) 

{Dt-aDl)h.f = 0, (84b) 

m f + agh = 0. (84c) 

We remark here that the bilinear equations (IM)l are presented in [l^ by using binary Bell polynomials 
method. To find multi-soliton solutions of Eqs. (IMll . we expand f,g and h as following 


/=l + e2/(2)+eV(4) + 
g = eg^^^ + + • • • , 

h = + e/i(3) + ... . 


(85) 

( 86 ) 

(87) 


Substituting the expansion (I55]) - (I571) into the bilinear Eq. (IM)) we find that there are only odd order terms of 
e in the first two equations while only even order terms appear in the third one. Comparing the coefficients 
at each order of e, we obtain one-soliton solution 


g = eexp(77), h = eexp(^), 


/ = 1 + 


—a 


Sa‘^{p + r)' 


■exp(ry + ^), 


where rj = px + qt + = rx + st + and q, s satisfy the dispersion relation 


s = ar 


q = -ap , 

Here p,r and are arbitrary constants. The two-soliton solution is presented as follows 

g = e(^exp{gi) + exp(r 72 )^ + e^(^ai 2 i exp{gi + g 2 + 6) + ai 22 exp(r 7 i + g 2 + 6)), 
h = e(^exp(^i) + exp(^2)) + (6121 exp(^i + 6 + 771) + 6122 exp(^i + ^2 + 772)), 

/ = 1 + cii exp( 77 i + ^1) + C12 exp( 77 i + ^2) + C21 exp{g 2 + ^1) + C22 exp{g 2 + ^2) 

+ e'‘ci 2 i 2 exp(? 7 i + g 2 + + ^ 2 ), 

where 


gi = piX + q^t + 77°, = TjX + Sjt + 

ft = -aph Sj = ar'j, 

_ _ Q!(pi -P2)^ __ ajpi -P2Y 

8a2(pi + ri)2(p2 + ri)2 ’ 8a2(pi + r2)2(p2 + 7 - 2 )^ ’ 

, __ ajri - r 2 Y , __ a(7-i - 7-2)^ 

So?{pi+riY{jpi+r 2 Y' 8a2(p2+ 7 'i) 2 (p 2+ 7’2)2’ 

c-. —__ 7 7 — 12 

8a2(p, + fe,)2’ 

^ _ _ c?{pi -P 2 )’^(ri - 7-2)^ _ 

64a4(pi + ri)(pi + r 2 )(j >2 + 7'i)(p2 + 7-2) ’ 

We can use the following compact expression for the two-soliton solution, 

g = exp{gi) + exp{g2) + a(l,2, l*)exp(77i -f 772 -f ^1) 

+ a(l, 2, 2*) exp(77i -b 772 -b 77*), 
h = exp(^i) -b exp(^2) + a(l, 1*, 2*) exp(77i -b 6 + 6) 

+ a(2,l*,2*)exp(772 +ei+6), 

/ = 1 -b a(l, 1*) exp(77i -b ^1) -b a(l, 2*) exp(?7i -b 6) 

-ba(2,l*)exp(772 -bCi) + a(2,2*)exp(772 -b6) 

+ a{l, 2,1*, 2*) exp(77i -b 772 -b Ci + 6), 


( 88 ) 

(89) 

(90) 

(91) 

(92) 

(93) 

(94) 

(95) 

(96) 

(97) 

(98) 

(99) 


( 100 ) 

( 101 ) 

( 102 ) 
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where the coefhcients are defined as 


a(.hj*) = - 
a{i,j) = - 


8a'^{pi + Tj) ’ 

8a'^{p^ -pjY 


i(* ,3 ) = - 


a 

8a? {vi — TjY 

5 

a 


^(^15^2;*'* 7^n) - 5 ^/) ■ 

l<k<l‘Cn 

In the same way, we give the exact A^-soliton solution of Eqs. ([51|) in the following form 

2N 


(e) 

/ = y] exp 
(o) 

g = exp 

L'^0,1 

(o) 

h = exp 

A=0,1 


N 2N 

j—1 j—N-\-l l<k<l 

N 2N 2N 

+ E '^3^3-N + E 

j—1 j—N-\-l l<k<l 

N 2N 2N 

?k^lAkl 

j—1 j—N-\-l l<k<l 


(103) 

(104) 

(105) 

(106) 


(107) 

(108) 

(109) 


where 


ip = PiX + q^t + pj qj =-apj j = l,2, • 

■■,N, 

(110) 

^j ~ I’jX Sj t , Sj = QTj , J = 1, 2, * • 

■ ,N, 

(111) 

( ^ \ 8a2(pfc -piY 

exp(Hfe;) = , k < 1 = 2,3,- ■ 

a 

■,N 

(112) 

O' 

exp{Ak,N+i) = „ n, , fc,/ = l,2,- 

8a^[Pk + n) 

■■,N, 

(113) 

(^ \ 8a2(rfc - riY 

exp{AN+k,N+i) = , k < 1 = 

2,3,--- ,N. 

(114) 


Here Pj , rj are both real parameters relating respectively to the amplitude and phase of the ith soliton. The 
summations X)jj=o n Sl=o i X)i=o i satisfy the condition (l26]) - (l28l) respectively. 


3.2 Integrable discretization of RNLS equation 

The differential-difference RNLS equation is obtained by discretizing the spacial part of the bilinear Eq. 


(115) 

(116) 

(117) 

(118) 

(119) 

( 120 ) 


^x9*f ^ g2(9n+lfn-l 2g„/„ -|-^n-i/n+i), 

where x = nS, n being integers and <5 a spatial-interval. We get 

DtPn • fn + ~^{9n-+lfn—l ~ ‘^9nfn + 9n—lfn+l) = 0, 
I^thn * fn (^n+1 fn —1 ‘^^nfn “b ^n—1 fn-t-1) — 0, 


8a^ 


{fn+lfn-1 - fn) + = 0. 


Let pn = e'^"' fn, hn = e"" fn- Then Eqs. (I116|) - (jll8|) transforms into the following nonlinear form 


'^n,t 


'^n.t CO 


( 1 - 




Un+Vn'^^gUn + l—Ur, ^Un-l—U 


(52 8a2 

n - ^ 

52 8 a 2 ‘ 


-be" 


‘)-2 






= 0 , 

= 0 . 


12 






















When we take the continuum limit (5 —?► 0, Eqs. (I119I) - (I120|) tend to Eqs. d7^ - (l80l) . In the sequel we 
use this discrete spacial step <5 = 1. Multiplying ()119|1 by 2(1 — + e""-!”""), and ()120|) by 

2(1 + a)(e“"+i““'* + adding and subtracting each other, respectively, yield 

2(1 + a) \un,t + a(e“"+i-“" + - 2)1 + e’'-!-''") 

+ 2(1 - a) \vn,t - a(e""+i-’'" + + 2)1 


2(1 - a) \un,t + a(e“"+i-“" + - 2)1 (e’'"+i-’'" + e’^-i-''") 

+ 2(1 + a) [vn,t - a(e""+^-"" + + 2)] 

_j_ ® g'Un + 'Un ^g'Un + l-Un _|_ g«„_ i^ ^gt;„ + i-n„ _|_ g!;„_i-i;„^ _ Q 


Adding and subtracting (I121|) - (jl22l) give us 


Un,t + a(e""+i-“’* + - 2) + e^'-i-"") 

Wn.t - a(e''”+i-’'’‘ + e’'"-!-"" + 2)1 (e“"+i-“" + e“"-i-“") = 0, 


4a 


un,t + 0(6“"+^-“" + - 2) + e’'-!-''-) 


- 4a vn,t - 0 ( 6 ’'"+!-"" + + 2) 

- ae“'‘+’'"(e“'‘+i"“" + e“"-i-“'‘)(e’^"+i“''” + e’'"-!"’'") = 0. 


From (I123D and (I124D . we get 


4(^ - a) u„,i + a(e“”+i-“" + - 2) (e^n+i-v,. g^n-i-«„) 

+ 4(a + i) \vn,t - + 2)1 (e“"+i-“" + 


ae 


y-n+Vn ^^Un + l—Un 




Setting 


TT ('^n \ '^n . Un \ -rr jj^ ( 

Un = exp ^-h la -^-J , V^ = U^= exp 


H” ' 

— la- 


( 121 ) 


( 122 ) 


(123) 


(124) 


(125) 


2 2/^ ^ ^V2 2/ 

or equivalently, 

Un = ln(C/„^“ ), = ln(t/„''“ 14''“ ), 

and substituting Un,Vn by Un, Vn into Eq. (I125p . we obtain a semi-discrete RNLS equation 

'a-iUnO / /TT^.,\^ / TA n \ AA 2 ^ ^ ^ ^ 2 / ^A , \ ^ 2 


(126) 


(127) 


a Un 
' a + i Un,t ' 


[-("f -(I(^)*r -1(^)*r -)i [i(^)^r -1^ 


a + ' 


Un 


/Un-l\ 2 a 
V Un ) 


/A„_1N■ 

V Un ) 


(i(^)^rH(%r)^r)(i(^)^rH(^)^i>»^ 


(128) 


In what follows we put a = 1 and a = —8 for the sake of simplicity. With these special parameters, the 
semi-discrete equations (I119I) - (I120I1 are simplihed as follows 


Un,t + (1 + e"“+’'“)(e“”+i-“" -b -2 = 0, 


'^n.t 


-(i + . 


Un+Vn^^^Vn + l—Vn 


_^g«„-l-«n)_^2 = 0. 


(129) 

(130) 


13 





























































































By the transformation e“" it„ and e’'" —>• Vm Eqs. (I129I) - (I130I) can be expressed in an alternative simpler 
form 


Un,t - ‘2Un + (1 + UnVn){Un+l + M„-l) = 0, 

Vn,t + 2n„ - (1 + UnVn)iVn+l + Vn-l) = 0. 


(131) 

(132) 


3.3 Soliton solutions of the integrable semi-discrete version of RNLS 

In this section, we construct the multi-soliton solutions of the semi-discrete RNLS equation (I128D and hence 
show its integrability. We rewrite Eqs. (I116I) - (I118|) with <5 = a = 1 in following compact form 


(^Dt - 2 + 2 cosh(i:)„^ g^,f„ = 0, 
(^Dt + 2-2 cosh(D„)^ hn»fn = 0 , 
2 sinh^ (^jfn»fn- gnhn = 0 . 


(133) 

(134) 

(135) 


Similar to the continuous case, we expand /mdri smd into series with respect to a small parameter e 
as follows 


/„ = 1 + + • • • + + '' • , 

+ . •. + + • • • , 

+ e3/,(3) ^ ^ + ... . 


We obtain the one-soliton solution 


2 / /3l + /?! 


= eexp( 77 i), = eexp( 77 (), /„ = 1 —csch (- ^ 


)exp(77i +-q[), 


where rji = ait + Pin -I- 71 , r][ = a'^t -I- P[n + 7 ^, and Pi, P'l satisfy the dispersion relation 


ai + 4sinh^(^) = 0 , 

a[ — 4sinh^(^) = 0 . 


Here 01,71 and oi, 7 ( are arbitrary parameters. The two-soliton solution is presented as follows 

gn = e[exp(77i) -b exp(772)] + e^[di exp(77i -b 772 -b ri'i) -b P 2 exp(77i -b 772 + ?72)]> 
hn = e[exp( 77 () -b exp(772)] -b e^[ci exp( 77 i -b 773 -b 771 ) -b C 2 exp( 77 i -b 773 -b 772 )], 

/„ = 1 -b ^csch^( ^^ ^^^ ) exp(77i -b ?7() -b ^csch^( ^^ ^^^ ) exp(77i -b 77 ^) 

+ ^csch2(^2-^)exp(772 -b77() -b icsch2(^^2-^) exp(772 -b 77 ^) 


+ e^x[exp(77i -b 772 -b 77 ^ -b 773 )], 


where 


Qi = 


^1 = 


_ g/32-l-/3i^2 

(e/3i+/31 _ l)2(e/32-H/31 _ 1)2^ 

^g/3i-l-/3i _ gPi+P' 2^2 


Q2 = 


(g/3i+/3; _ l)2(g/3i+/3' _ ;^)2 ’ I'g/Ss-S/SJ _ ;^)2^g^2-S/3' _ ly ’ 

_ g^2^2g/3i-)-/32+/3j+/32 

(e/32+^1 _ l)2('g^2+/3^ _ l)2(e^i+/31 _ l)2(g/3i-h^^ _ * 


^g/3i+/32 _ g/32+/32^2 

(e/3i+^^ _ l)2(e/32+/3^ _ 1)2 ^ 
^g^2+/3i _ g/32+/32^2 


X = 


(136) 

(137) 

(138) 


(139) 

(140) 

(141) 


(142) 

(143) 


(144) 


(145) 
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Fig. 3: The resonant interaction of two solitons of discrete RNLS. One soliton splits into two solitons, or, 
two solitons fuse into one. 


The coefficient x plays a role of classification of soliton interactions. When x is a finite number and not 
equal to zero, the regular soliton interaction exists. The resonant interactions of two solitons occur in Eq. 
(I128|) in the case of x —>■ oo or x —>■ 0. Our discussion is focused on the case of x = 0 since the other case 
can be analyzed in a similar way. Parametric conditions is given as follows 

= q. (146) 

Choosing different factors of (11461) as zero, we get two types of soliton resonances, i.e., the fission and fusion. 
FigEl (a) describes one soliton breaks into two solitons with the parameters selected as /3i = 0.9 ,/32 = 
1.8,/3( = 0.3 ,/32 = 0.3. Figl3](b) shows two solitons fuse into one soliton with the parameters /3i = 0.8 ,/32 = 
0.8,/?( = 1.5,= 0.9. 

Similar as the two-soliton solution (I100D - (I102I) for the continuous equation, we can use the following 


compact expression for the above two-soliton solution, 

/„ = 1 -b a(l, 1 *) exp{r]i + r][) + a(l, 2 *) exp{r]i + ij'^) 

+ a( 2 , 1 *) exp(?72 -b r][) + a( 2 , 2 *) exp(?72 -b rj'^) 

+ a(l, 2,1*, 2*) exp(? 7 i + m + Vi + V 2 ), (147) 

g„ = exp(? 7 i) -b exp(? 72 ) -b a(l, 2 , 1 *) exp(? 7 i + r ]2 + r][) 

+ a(l, 2 , 2 *) exp( 77 i -b 772 + V2), 
hn = exp{r][) + exp{r] 2 ) + a(l, 1*, 2 *) exp( 77 i + r][ + 

+ a(2,r,2*)exp(772+77l+772), (148) 

where the coefficients are defined by 

a{i,j) = 4sinh^( ^' ^ (149) 

a{hj*) = ^csch^(^^-^), (150) 

a{i*,j*) = 4sinh^( ^' ^ (151) 
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The exact A^-soliton solution of eqs. (I133I) - (I135I) is presented in the form 


where 


(e) 

N 

2N 


2N 



/n = ^ exp 



+ 


(152) 


_i=i 

j=N+l 


l<k<l 



(o) 

N 

2N 


2N 



= X] 

+ 



y ^ vkviAki 

(153) 


J = 1 

i=N+l 

l<k<l 



(o) 

N 

2N 


2N 



hn = y] exp 

y . + 

Y ^3'n'o-N 

+ 

y ^ XkXiAki 

(154) 

A^0,1 

1=1 

j=N+l 


l<k<l 



r/j = ajt + PjTi 

+ 71, 

= _4sinh2^, 
2 


J = l,2,- 

■,N, (155) 

Tj'j = a'jt + PjTi 

+ 7j, 

2/3^- 

= 4sinh 

2 

j 

= 1,2,... 

,N, 

(156) 

exp(Afcz) = 4sinh^(^^ ^ 

-), k<l-2 

3 

••• ,iV, 


(157) 

exp(Afc,Ar+i) = 

^CSCh^i^'' 

4 

k.i- 

1, 

2,--- ,N, 


(158) 

exp(AAr+fe,7V+i) - 4smh ( 


= 2,3,--- 

,iV. 

(159) 


The summations EL°io,i ^nd 




4 Conclusions 

To summarize, we presented here one semi-discrete integrable version for the (2-1-1 )-dimensional modified 
HF equation and one semi-discrete version for the resonant NLS equation. TV— soliton solution to the both 
discrete systems are given in standard form. Interaction properties of two-soliton solutions are analyzed. 
The regular, intermediate-state and resonant soliton interactions of the semi-discrete systems are discussed. 
Those analysis on the resonant interactions might have the applications in optical communication systems. 
Fully discrete integrable version for the two equations is under consideration. 
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